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Abstract: The phase diversity technique is a useful tool to measure and 
pre-compensate for quasi-static aberrations, in particular non-common path 
aberrations, in an adaptive optics corrected imaging system. In this paper, 
we propose and validate by simulations an extension of the phase diversity 
technique that uses long exposure adaptive optics corrected images for sens- 
ing quasi-static aberrations during the scientific observation, in particular for 
high-contrast imaging. The principle of the method is that, for a sufficiently 
long exposure time, the residual turbulence is averaged into a convolutive 
component of the image and that phase diversity estimates the sole static 
aberrations of interest. 

The advantages of such a procedure, compared to the processing of short- 
exposure image pairs, are that the separation between static aberrations and 
turbulence-induced ones is performed by the long-exposure itself and not 
numerically, that only one image pair must be processed, that the estimation 
benefits from the high SNR of long-exposure images, and that only the static 
aberrations of interest are to be estimated. Long-exposure phase diversity 
can also be used as a phasing sensor for a segmented aperture telescope. 
Thus, it may be particularly useful for future planet finder projects such as 
EPICS on the European ELT. 
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1. Introduction 

Calibrating the quasi-static aberrations of a ground-based Adaptive Optics (AO) [1] corrected 
imaging system is an important issue, especially for extreme AO high contrast instruments such 
as the proposed planet finder instruments for the ESO and Gemini 8-meter telescopes. 

Measuring these aberrations off-line, i.e., during a day-time calibration on an internal refer- 



ence source, has been successfully applied to existing systems such as Naos [2, 3] or Keck, 
and recently refined for the SPHERE project in order to achieve a nanometric accuracy [4]. 
The main limitations of such a procedure directly stem from its off-line nature: the aberrations 
located before the internal reference source are not sensed, and the aberrations may evolve be- 
tween the day-time calibration and the night-time observation. These two limitations can be 
circumvented by an appropriate calibration performed on-line, i.e., during night-time observa- 
tions, as proposed in the following. 

Two problems must be addressed in order to calibrate on-line the quasi-static aberrations 
of the optical system made of the telescope, its AO system and the camera: the first one is 
to distinguish between the turbulence-induced component of the wavefront and the static one, 
which is the only wavefront of interest for the problem at hand, and the second one is to sense 
all the aberrations from the (potentially segmented) primary down to the focal plane of the 
camera. 

A wave-front sensor (WFS) is able to measure the aberrations seen by the telescope on-line, 
but these consist of the sum of a turbulence-induced component, which is partially compen- 
sated for by the AO, and a static component. Because the turbulence, whether corrected or not, 
evolves quickly, the WFS measurements are generally performed with integration times that 
freeze the turbulence evolution, typically a few milli-seconds. 

There is today a large number of WFSs, which are thoroughly reviewed in Ref. [5] and can be 
classified into two families: pupil-plane sensors, such as the Hartmann-Shack and the curvature 
sensors, and focal-plane sensors. A focal-plane WFS is the only way to be sensitive to all aber- 
rations down to the focal plane, and in particular to the so-called non-common path aberrations 
of an AO system, which motivates our choice for a focal-plane WFS in the following. Estimat- 
ing aberrations from a single focal-plane image of a point source is a difficult problem known 
as phase retrieval. Phase-retrieval has two major limitations. Firstly, it only works with a point 
source. Secondly, there is generally a sign ambiguity in the recovered phase, i.e., the solution is 
not unique. Gonsalves [6] showed that by using a second image with an additional known phase 
variation with respect to the first image such as defocus, it is possible to estimate the unknown 
phase even when the object is extended and unknown. The presence of this second image addi- 
tionally removes the above-mentioned sign ambiguity of the solution. This technique, referred 
to as phase diversity, has been significantly developed in the past twenty years, both for wave- 
front sensing and for imaging; see for instance Refs. [7, 8, 9, 10, 11, 12, 3], and Ref. [13] for a 
review. 

The estimation of static aberrations from a series of short-exposure phase-diversity data has 
been performed using a series of image pairs of an astronomical object [14, 15]. Lee et al. [16] 
have performed such a calibration of static aberrations with a series of images instead of im- 
age pairs, and an original diversity: no additional defocused image was used, and successive 
changes to the adaptive optics introduced the required diversity. 

In both approaches, the static aberrations are obtained as an empirical average of the phase 
estimates corresponding to each short-exposure data. This is notably suboptimal for at least 
three reasons: 

• the images correspond to short integration times, and are consequently noisier than the 
corresponding long-exposure image pair, so each phase estimate suffers from this noise; 

• the computational cost is high because many short-exposure images must be processed 
in order to estimate the sought static aberrations; 

• the phase estimation accuracy may be penalized because the estimation must be per- 
formed on a number of phase parameters that is large enough to describe the short- 
exposure phase, whereas only a smaller number of these parameters may be of interest. 



if the sought static aberrations are of lower order than turbulence-induced ones; 



In this paper, we propose and validate an extension of the phase diversity technique that 
uses long-exposure AO-corrected images for sensing quasi-static aberrations. This way, (1) 
the separation between quasi-static aberrations and turbulence-induced ones is performed by 
the long-exposure itself and not numerically, (2) only one image pair must be processed, (3) 
the estimation benefits from the high SNR of long-exposure images, and (4) only the static 
aberrations of interest are to be estimated. 

2. Principle of long-exposure phase diversity 

We consider a ground-based telescope observing Space through the turbulent atmosphere. The 
long-exposure optical transfer function (OTF) of the atmosphere-ninstrument system is the prod- 
uct of the OTF of the sole instrument /is, called static OTF in the following, by the atmospheric 
transfer function (ATF) [17]: 

{ho)=kk. (1) 

The static OTF is a function of the unknown static aberrations, which are coded in the phase 
function (p in the aperture; let P be the indicator function of the aperture, i.e., 1 in the aperture 
and outside, is given by: 

h,{(p)^Pe"f®Pe''<' (2) 

where ® denotes auto-correlation. The phase function <p{u,v) is expanded on a basis {bk}, 
which is typically either Zernike polynomials or the pixel indicator functions in the aperture: 
^(m,v) = Y.k'Pkbkiujv), where the summation is, in practice, limited to the number of coeffi- 
cients considered sufficient to correctly describe the static aberrations to be estimated. We shall 
denote by (j) the vector concatenating the set of unknown aberration coefficients 0^^. 

Assuming phase perturbations with Gaussian statistics, the ATF at any spatial frequency f is 
given by [17]: 

^^(f) =e-i^'/>(^f) (3) 

where A is the imaging wavelength and is the phase structure function. If the turbulence 
is partially compensated by an AO system. Equations (1) and (3) remain valid [1], although 
slightly approximate because the residual phase after AO correction is not stationary. 

With these equations in hand, and assuming that the image is not larger than the isoplanatic 
patch, we can now model the long-exposure image. The image is recorded by a detector such as 
a CCD camera, which integrates the flux on a grid of pixels. This can be conveniently modeled 
as the convolution by a detector PSF h^, assumed to be known in the sequel, followed by a 
sampling operation. Using Eq. (1), the global long-exposure PSF of the instrument is thus: 

hie = hi ★ {ho) = /id * K * h^, (4) 

where /zs is the PSF due to static aberrations, given by the inverse Fourier transform of Eq (2), 
/za is the atmospheric PSF given by the inverse Fourier transform of Eq. (3), and ★ denotes 
convolution. 

Due to the inevitable sampling and noise of the detection processes, the image if recorded 
in the focal plane is the noisy sampled convolution of the long-exposure point-spread function 
(PSF) /lie with the observed object o. This model is generally approximated by a noisy discrete 
convolution with the sampled version o of the object o: 



if~hif.^o+n, (5) 

where hi^, is the sampled version of /iie, and « is a corruptive noise process. If the noise is not 
additive and independent from the noiseless image, for instance if it is predominantly photon 
noise, then Eq. (5) should read i —h\^-koon, with the symbol o representing a pixel-by-pixel 
operation [18]. For legibility we shall keep the additive notation. 

Combining Eq. (5) with the discrete counterpart of Eq. (4) yields the following discrete image 
model for the focused and defocused images respectively: 

If = hi-kh^^\{hai<o)+n (6) 
Id = hii.h^^^'^'K{Ki.o)+n', (7) 

where /id, and are the sampled versions of /id, /Js and h^, (pd is the known additional phase 
introduced in image I'd, and the superscripts on are reminders of the aberrations that enter the 
static PSF of each image. 
Let o' be the convolution of the atmospheric PSF with the observed object: 

o' = h^i^o (8) 

The phase-diversity data model of Eqs. (6) and (7) is strictly identical to the one that would 
be obtained by imaging the pseudo-object o' in the absence of turbulence and with the same 
static aberrations. Thus, all the phase estimation methods developed for short-exposure images, 
in which the OTF of the system is completely described by a phase function, can be applied 
here to estimate the sole static aberrations. 

Additionally, it is well-known that, for a given noise level, the estimation quality of the 
aberrations in phase diversity depends on the spectral content of the observed scene. Thus, in 
the method proposed here, the estimation quality of the aberrations will depend both on the 
spectral content of the observed object and on the ATE, i.e., on the turbulence correction quality 
provided by the AO. 

The appropriate implementation of this long-exposure phase diversity technique depends on 
the type of instrument. For non coronagraphic instruments, one should use images provided by 
the science sensor; the defocused image can be either obtained simultaneously with the focused 
science image by means, e.g., of a beamsplitter, or alternately. In the latter case, the deformable 
mirror itself can be used to provide the defocus [2, 3]. In both cases the fraction of the incoming 
flux allotted to the defocused image may be notably less than 50%, in order to maximize the 
flux on the scientific data. Indeed, if the quasi-static aberrations are measured at intervals of, 
e.g., a half-hour, defocused images must only be available with that kind of rate. 

For coronagraphic systems, for which the aberrations to be minimized are the ones located 
before the coronagraph, one may use a beam-splitter and an auxiliary image sensor located 
just before the coronagraph. For the SPHERE instrument such a sensor actually already exists 
in the design for centering the star image on the coronagraph: it is the so-called differential 
tip-tilt sensor. This sensor could easily be adapted and used for the long-exposure phase diver- 
sity measurements. As in the non coronagraphic case, the defocused image may be obtained 
simultaneously with the focused image by use of a beamsplitter, or alternately by means of a 
longitudinal displacement of the sensor by a few millimeters. 

3. Chosen phase estimation method 

Among all the possible estimation methods (see, e.g., Ref. [13] for a review) in this paper we 
choose, for simplicity, the conventional least-squares joint estimation of the phase (/) and the 
object o', with a regularization on both quantities: (d', 0) = argmin7(o', 0) with 



Jio\^)^-^\\if-hi*hi^Ko'\\ + -^\\i,-h,*hi^+^''^^o'\\+R,,{^^^ (9) 

where and CJy are the noise variances of the two images, estimated beforehand. The object 
regularization is chosen as a quadratic function, so that the whole criterion J is quadratic with 
respect to o' and thus has an unique, closed-form solution o'{(p) for a given phase 0. This 
allows one to replace the optimization of J{o', (p) with that of criterion J'{(p) = J{o'{(j)), (p), as 
commonly done in the unregularized case [6, 8]. 

Following the findings of Blanc [19, 12], we under-regularize the object in order to best 
estimate the phase. This strategy is supported by the fact that it yields a phase estimation with 
satisfactory asymptotic properties, as shown in Ref. [20], and that these properties hold even if 
the noise is not Gaussian. 

Concerning the phase, we choose the basis of the pixel indicator functions rather than, e.g., 
a truncated basis of Zernike polynomials, in order to model and reconstruct phases with a high 
spatial frequency content. Because of the potentially large number of phase unknowns we are 
lead to regularize the phase estimation. To this aim, we use a functional proposed specifically 
for such a phase basis in Refs. [19, chap. 7] and [13, Sect. 8], which is recalled below: 



J('Pl.m-l-<l>lji,) _ gj{<l>l.m-<l>lji,+ l) 

where the summation is done on all the pixels within the pupil {S is the pupil support). Further- 
more, we impose a strict support constraint i.e., all terms that contain, at least, a pixel out 
of the pupil support, are suppressed. 

This regularization function has been constructed in such a way that it is insensitive, as the 
data is, to a global piston, to tip-tilt, and to any In variation of the phase on any pixel. This 
way, no local minimum is introduced by the regularization into the minimized criterion. 

4. Validation by simulations 

We shall now validate the proposed method by simulations. We shall essentially study the in- 
fluence of the exposure time. Indeed, the main specific assumption of the proposed method lies 
in Eq. (1), because the factorization of the OTF in a static OTF and an ATF is strictly valid only 
if the turbulence is perfectly averaged i.e., for an infinite exposure time. Note that by exposure 
time we mean the (finite) number of independent turbulence realizations, not the noise level. 
In all the simulations presented here, we have considered noiseless images. We have checked 
that the behavior of the phase estimation in the proposed method with respect to the noise level 
is not specific and is the same as conventional phase diversity with short-exposure images: the 
estimation error is usually proportional to the average standard deviation of the noise in the 
images [21, 12]. Then we shall briefly study the influence of the AO correction quahty. 

4.1. Conditions of simulation 

We consider here a point-like source, observed with an 8 m ground-based telescope, equipped 
with AO. The simulations take into account both the AO-corrected turbulence and static aber- 
rations. The baseline adaptive optics system considered is the high-performance AO system 
Saxo [22] of the SPHERE [23] instrument. We use a Fourier-based simulation method that 
describes the AO via the spatial power spectrum of the residual phase [24] and is presented in 
Ref. [25]. The simulation takes the following reaUstic set of parameters: a 41 x 41 subaperture 



f,j[<i>i-i., 



„j{^Lm-<l>l+l.m) 



Shack-Hartmann, a 1.2 kHz sampling frequency, a guide star of magnitude 8, and a Paranal- 
like turbulence profile, with a seeing of 0.8 arcsec at 0.5 [im. Static aberrations are randomly 
generated according to a /^^ spectrum, / being the spatial frequency in the pupil, with a total 
wavefront error of 0.23 radian RMS at 1.6|jm, i.e., 60 nm RMS. 

For some of the simulations, the number of actuators on a pupil diameter will be decreased 
from 40 X 40 to 14x14 and 7 x 7 in order to study the influence of the AO correction quality. 
The above simulation conditions lead to a phase variance of the residual turbulence which is 
respectively 1, 7 and 21 time(s) the variance of the static aberrations. 

Two simultaneous long-exposure images are simulated, with a phase diversity between these 
images consisting of a 1.814 radian RMS defocus. The corresponding defocus distance is pro- 
portional to the square of the f-number of the system and, at 1.6[im, is 2.9mm for an //15 
system such as Naos and 2cm for an //40 system such as SPHERE. These images are simu- 
lated following two schemes; 

• Finite exposure time images are made of the summation of short exposures. The long- 
exposure OTF is then the sum of short-exposure OTFs, each of which being computed 
through Eq. (2), with a phase (p composed of the sum of the static aberrations and of the 
instantaneous AO-corrected turbulent wavefront. The turbulent wavefront is randomly 
generated from a PSD that takes into account both the turbulence profile and the AO 
correction [25]. 

• Infinite exposure time images are not computed as an empirical average. Instead, the 
residual phase structure function is computed from the PSD of the AO-corrected 
turbulent wavefront, then the ATF is computed via Eq. (3), and the images if and are 
computed according to Equations (6) and (7). 

4.2. Influence of exposure time 

We first consider the case of a high performance, SAXO-like correction. Figure 1 shows the 
evolution of the reconstruction error with the number of exposures used for the simulation. The 
first points (from 10 to 1000 exposures) are simulated with a finite exposure time, whereas the 
last point (noted infinity on the X-axis) is simulated with an infinite exposure time. Because 
the correlation time of corrected turbulence is typically 10 to 100 milliseconds depending on 
turbulence parameters and on the AO correction quality, the simulations with 1000 exposures 
correspond to an integration time between 10 and 100 seconds. 

The reconstruction error decreases with the number of exposures, down to very weak values 
(less than 0.01 radian) for an infinite exposure time. As the number of exposures increases, 
Eq. (1) becomes more valid, turbulence residuals become better fitted by a modification of 
the estimated object only (into an object o' ~h^*o, cf. Eq. 8), and the estimated phase is 
eventually only the set of static aberrations. Incidentally, we have checked that for a single 
turbulent exposure, the phase estimated is the sum of the turbulent wavefront and the static 
aberrations. 

This simulation study shows that the estimation of static aberrations from a single pair of 
turbulence-degraded images is possible. The quality of aberration reconstruction is directly 
linked to the convergence of the images towards long-time exposures. For the AO system and 
the level of static aberrations considered here, an integration time corresponding to a thou- 
sand independent turbulence realizations yields a phase estimation error of about 0.012 radian, 
close to that obtained with an infinite exposure time. At 1.6fim, this number translates into a 
3nm RMS optical path difference. This precision is compatible with the 5 — lOnm RMS static 
aberration residual that is needed for the detection of warm Jupiters on an 8-meter telescope. 
Incidentally, we see that with less than 1000 exposures the required precision would not be 
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Fig. 1. Evolution of the reconstruction error with the exposure time (in number of inde- 
pendent turbulence realizations). Static aberrations are randomly generated according to 
conditions described in Subsection 4.1. AO correction is assumed to be performed by a 
SAXO-like system. 



obtained for such a mission. 

4.3. Influence of AO correction quality 

Figure 2 shows the evolution of the estimation error, for different correction qualities, obtained 
here simply by varying the number of actuators of the AO system, all other parameters being 
equal. One can see that as the correction quality degrades, the turbulence residuals are more im- 
portant and thus the estimation of the static aberrations needs more exposures for the same error 
level. The estimation errors for an infinite exposure time are equivalent for all three correction 
qualities. 

In the case of a Naos correction with 14 x 14 actuators and 1000 exposures, the phase es- 
timation error is about 0.13 radian, which at 2.2|jm translates into a 45nm RMS optical path 
difference. This is almost three times smaller than the residual static aberrations of 120nm RMS 
measured on Naos-Conica after off-line phase diversity measurement and correction [3]. The 
on-line long-exposure phase diversity technique could thus be an attractive way to calibrate 
quasi-static aberrations on non-extreme AO systems too. 

We now detail the spectral analysis of the estimated aberrations. In an AO system, the number 
of actuators determines the highest spatial frequency of the turbulence to be corrected. This 
parameter therefore directly impacts on the spectral content of the turbulence residuals. On 
the following figures we plot the circularly averaged spectra of the estimated aberations with 
respect to the spatial frequency, for different AO correction levels. 

Figure 3 shows the error spatial spectrum in the case of a Naos correction (7x7 actuators), 
for different exposure time. One can see that the spectrum of estimated aberrations get closer to 
that of the true ones as the number of exposures increases, and that the convergence is slower 
for the uncorrected frequencies of the turbulence (above 3.5 cycles/pupil for this 7x7 actuator 
system). 

In the case of a higher correction (respectively 14 x 14 and 40x40 actuators, on Figures 4 
and 5), the turbulence residuals are reduced and the convergence of the spectrum to true 
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Fig. 2. Evolution of the reconstruction error with the exposure time, for several levels of 
AO connection. Static aberrations are randomly generated according to conditions described 
in Subsection 4. 1 . 



spectrum is consequently faster. Moreover, as in the Naos-7 case, convergence is notably 
faster in the corrected part of the aberration spectrum, which correspond to a limit of 7 and 
20 cycles/pupil respectively. 
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Fig. 3. Spatial spectrum of the estimated aberrations, in the case of a Naos-7 correction. 
The vertical line represents the maximum spatial frequency that is corrected by the AO 
system. 
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Fig. 4. Spatial spectrum of tiie estimated aberrations, in thie case of a Naos-14 correction. 
The vertical line represents the maximum spatial frequency that is corrected by the AO 
system. 
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Fig. 5. Spatial spectrum of the estimated aberrations, in the case of a Saxo/Sphere cor- 
rection. The vertical line represents the maximum spatial frequency that is corrected by the 
AO system. 



5. Conclusion and perspectives 

The phase diversity technique is a powerful tool to measure and pre-compensate for quasi-static 
aberrations, in particular non-common path aberrations, in an AO-corrected imaging system. 
So far it has, to the best of our knowledge, either been used off-line on an internal calibration 



source, or on-line from short-exposure turbulence-degraded data. Each approach has its own 
limitations, discussed in this paper We have proposed and validated by simulations an extension 
of the phase diversity technique that uses long exposure AO corrected images for sensing quasi- 
static aberrations. 

The principle of the method is that, for a sufficiently long exposure time, the residual turbu- 
lence is averaged into a convolutive component of the image and that phase diversity estimates 
the sole static aberrations of interest. 

Technically, the advantages of such a procedure, compared to the processing of short- 
exposure image pairs, are that the separation between static aberrations and turbulence-induced 
ones is performed by the long-exposure itself and not numerically, that only one image pair 
must be processed, that the estimation benefits from the high SNR of long-exposure images, 
and that only the static aberrations of interest are to be estimated. Compared to pupil-plane 
wavefront sensing techniques, phase diversity has the advantages that it senses all aberrations 
down to the focal plane and that the hardware is extremely simple. 

From a system point of view, on-line long-exposure phase diversity opens a new area of 
applications, in particular it will allow one to correct in real time (meaning during the scientific 
exposure) for any evolution of instrumental defects. This may be considered for improvements 
to the Sphere instrument and should significantly improve the overall system detectivity. 

This technique can also be used as a phasing sensor for a segmented aperture telescope. 
Indeed, phase diversity can be applied to the peculiar aberrations constituted by the differential 
tip-tilts and pistons of such a telescope. 

Thus, long-exposure phase diversity may be particularly useful for the future Planet Finder 
project on the E-ELT called EPICS [26]. Indeed, on the one hand, for this project, the detectiv- 
ity requirements are by far more stringent than for SPHERE and the on-line correction of non- 
common path aberrations is mandatory. And on the other hand, without any opto-mechanical 
modification of the sensor, long-exposure phase diversity should enable very accurate measure- 
ments of the segments' phasing. 

Lastly, for some applications it may be useful to estimate not only the static aberrations but 
also the ATE, which is part of the estimated quantities. The method proposed herein can be 
refined to estimate the ATE precisely using (a) a parametrization of the ATE through the phase 
structure function and (b) a specific space-varying regularization criterion for the ATE along 
the lines of [27, 28]. 

Short-term perspectives include a more thorough study of the performance of the long- 
exposure phase diversity technique, coupled with a global system analysis. 



